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Abstract 

We study charge k SU{2) BPS monopoles which are symmetric under the cyclic 
group of order k. Approximate twistor data (spectral curves and Nahm data) is 
constructed using a new technique based upon a Painleve analysis of Nahm's equation 
around a pole. With this data both analytical and numerical approximate ADHMN 
constructions are performed to study the zeros of the Higgs field and the monopole 
energy densities. The results describe, via the moduli space approximation, a novel 
type of low energy k monopole scattering. 
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1 Introduction 



By studying rational maps it has been shown || that the imposition of cyclic Ck symmetry 
upon strongly centred charge k SU{2) BPS monopoles selects out from the general k- 
monopole moduli space a set of k 2-dimensional totally geodesic submanifolds. Each 
of these 2-dimensional submanifolds is a surface of revolution, hence by application of 
a reflection symmetry a number of geodesies in Aik are obtained. In the moduli space 
approximation |§ these geodesies may be interpreted in terms of the motion of k slowly 
moving monopoles. It is therefore of considerable interest to know more about these one- 
parameter families of monopoles. The rational map approach is limited for this purpose, 
and the only additional information it supplies concerns the asymptotic in and out states 
for the scattering process. 

Other twistor approaches yield more information about the monopole, but are more 
difficult to apply than the rational map method. In || only the form of the corresponding 
spectral curves was found. In this paper we do a little better, by constructing approximate 
spectral curves, which of course have this required form. Moreover, we construct approxi- 
mate Nahm data and implement an ADHMN construction to compute the Higgs field and 
energy density of these monopoles. 

The technique used is to make use of the integrability of Nahm's equation, not to 
attempt an explicit solution, but rather to study the nature of a solution around a given 
pole, with the aid of a little Painleve analysis. The method applies equally well to all 
values of k, but only the cases k = 3 and k = 4 will be discussed in detail. 

2 Monopoles and ADHMN 

SU (2) BPS monopoles are topological soliton solutions of a Yang-Mills-Higgs gauge theory 
with no Higgs self-coupling. Boundary conditions imply that the Higgs field at infinity 
defines a map between two-spheres. This map has an integer valued winding number k, 
which we identity as the magnetic charge or number of monopoles. In this paper we shall be 
concerned with strongly centred monopoles, which basically means that we fix the centre of 
mass of the monopole configuration to be the origin, and set a total phase to be unity. For 
more precise details of strong centring and other background information on monopoles we 
refer the reader to @, 0]- 

As mentioned in the introduction, requiring invariance of a fc-monopole under cyclic 
Cfc symmetry, and an additional reflection symmetry, leads to a number of geodesies T} k 
(we follow the notation of 0) in the k- monopole moduli space M.k- Essentially there are 
(2fc + 3 + (— l) fc )/4 different types of these cyclic geodesies, corresponding to / = 0, 1, ..k/2 if 
k is even and I = 0,1, ..(k — l)/2 if k is odd. Physically, for I ^ 0, the associated monopole 
scatterings are distinguished by having the out state (or in state by time reversal) consisting 
of two clusters of monopoles with charges k — I and /. This explains why we do not allow 
I > k/2, since this is basically the same scattering event as one of the geodesies with I < k/2. 
If I = then the monopoles remain in a plane and scatter instantaneously through the 
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axisymmetric fc- monopole and emerge with a ir/k rotation. This kind of it/ k scattering is 
essentially a two-dimensional process and has been extensively studied in planar systems 
||. In this paper we shall be concerned with the more exotic scatterings with 1^0. It 
should also be pointed out that the case k = 2 is special, in that the geodesies S° and 
E2 are isomorphic, so that there is only one type of scattering; the famous right-angle 
scattering found by Atiyah and Hitchin ||. 

In this paper we shall only deal with the cases k = 3 and k = 4. The method applies 
equally well to all values of k, but all the important features are captured by these two 
examples. From the above we see that for k = 3 there is only one interesting geodesic E3. 
It describes three individual monopoles which scatter and emerge as a single monopole and 
a 2-monopole cluster. It also contains the tetrahedral 3-monopole || |6| as an instantaneous 
configuration. For k = 4 there are two interesting geodesies £4 and T*\. The first describes a 
scattering in which the monopoles emerge as two 2-monopole clusters, and also includes the 
cubic 4-monopole |6|. The second describes a scattering which results in the formation 
of a single monopole and a 3-monopole cluster; no monopoles with the symmetries of a 
Platonic solid are contained in this geodesic. 

There are several twistor techniques which are applicable to monopoles, but here our 
main tool will be the ADHMN |Hj [| construction. This is an equivalence between k- 
monopoles and Nahm data (Ti,T2,T3), which are three k x k matrices which depend on a 
real parameter s G [0,2] and satisfy the following; 

(i) Nahm's equation 



(ii) Tj(s) is regular for s G (0, 2) and has simple poles at s = and s = 2, 

(iii) the matrix residues of (T 1 ,T 2 ,T 3 ) at each pole form the irreducible fc-dimensional 
representation of SU(2), 

(iv) Us) = -T}(s), 

(v) T 4 ( a ) = 2?(2-*). 

Equation (i) is equivalent to a Lax pair and hence there is an associated algebraic curve, 
which is in fact the spectral curve 0j]. Explicitly, the spectral curve may be read off from 
the Nahm data as the equation 

det(77 + (7\ + iT 2 ) - 2tT 3 ( + (T t - iT 2 )C 2 ) = 0. (2.2) 

It is useful to construct the spectral curve since it not only gives a convenient representation 
of the monopole, but also furnishes the constants for Nahm's equation. 
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The procedure by which the Higgs field (and gauge potential) can be reconstructed 
from the Nahm data is outlined in section 5, where our analytical approximation to this 
procedure is also discussed. 

In the next section we give the general form for cyclic /c-monopole Nahm data and 
introduce the relevant equations for the case k = 3. 



3 Cyclic Nahm data 



As explained in |L3[ the form of the Nahm data for C& symmetric fc-monopoles may be 
obtained as a linear sum of generators of the affine Lie algebra Explicitly, consider 



the Lie algebra A 



0,..,k — 1 the generators of the extended Cartan 



subalgebra and E±i the generators corresponding to the simple roots at 



l,..,k-l, 



plus the lowest root ao = — Y^=\ oij. In the Chevalley basis these satisfy 



[Hi, E±j] 
[Ei, E_j] 



±CijE±j 
5ijHj 



where are the elements of the k x k extended Cartan matrix given by 



2(a i ,a j ) 



1,3 



0,..,k-l. 



(3.1) 
(3.2) 



(3.3) 



Cyclic Nahm data may be expressed in terms of a linear sum of generators as 

fc-i 



fc-i 

Ti -iT 2 = J2Qj E +ji 

j=0 



3=0 



(3.4) 



with real function coefficients Pj,Qj. 

For the case k = 3 a change of notation yields the Nahm data for C3 symmetric 3- 
monopoles to have the form 






-h 


h 




' 


h 


h ' 


h 





-h 




h 





h 


-h 


h 







h 


h 






h 
h 
-(h + fc) 



with corresponding spectral curve 

?? 3 + ar/C 2 + 7C 3 + ^(C 6 + i) 



0. 



(3.5) 



(3.6) 
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The constants a, j3, 7 can all be taken to be real, by imposition of a reflection symmetry, 
and are given in terms of the /j's as 



<*= fl + fl + fl-^fl + fl + Uh) 

P = -A/2/3 

7= 2/ 1 2 /4 + 2/ 2 2 /5 - 2(/ 3 2 + 4/ 4 / 5 )(/ 4 + / 5 ). (3.7) 
For this data Nahm's equation becomes 

fi= (2/ 5 + / 4 )/i 
k= -(2/ 4 + / 5 )/ 2 

/3 = (/4 — /5J/3 

A= (/ 3 2 -/ 2 2 )/2 

h= (/! 2 -/ 3 2 )/2. (3.8) 

The first calculation we require is to see how the tetrahedral 3-monopole sits inside this 
Nahm data. The spectral curve and Nahm data of the tetrahedral 3-monopole have been 
computed H for the monopole in a different orientation and it is possible to derive the 
answer we require by rotation of this known data. However, this is a non-trivial exercise 
which in fact requires more work than simply solving the equations again. Thus in the 
remainder of this section we compute the Nahm data and spectral curve of the tetrahedral 
3-monopole in the orientation in which it has C3 symmetry around the 23 axis. 

For tetrahedral symmetry we must have that a = 0, which is achieved by the following 
choice 

/ r = 2/ 4 2 , fi = 2/|, fl = 2(/ 4 + / 5 ) 2 . (3.9) 
The remaining constants then simplify to 

7 = -20/ 4 / 5 (/4 + / 5 ), and = (3.10) 

Using the above constants we can solve for / 5 in terms of / 4 as 



2/ 5 = -/ 4 ± V/ 4 2 - 7/(5/4) (3.11) 

and substituting this into the equation (|3.8|) for / 4 yields 

/ 4 2 = / 4 (/ 4 3 - 7/5). (3.12) 

Introducing the variable 2w = — 1//4 this equation becomes the standard form elliptic 
equation 

Aw 2 = cw 3 + 1, = 87/5. (3.13) 

Set 

w = p(t)/(4 K ) (3.14) 
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with t = ks + B then, choosing An = c 1 / 3 , p is the Weierstrass elliptic function satisfying 

p' 2 = V + 4 (3.15) 

where prime denotes differentiation with respect to the argument. The real period of this 
elliptic function is 

2wi = r(l/6)r(l/3)/(2^r). (3.16) 

For the correct boundary conditions take the elliptic function to go between zeros at t — 
2uoi/3 and t = 4a>i/3. This requires that k = cji/3 and B = 2ui/3. Now 

-0 = -^ = -£==^ = r(l/6) 3 r(l/3) 3 /( 7 r 3/2 3 3 8V2). (3.17) 

Thus we have arrived at the spectral curve 

rf - t(( 6 + i5V2( 3 + l)r(l/6) 3 r(l/3) 3 /(7r 3 / 2 3 3 8v / 2) = 0. (3.18) 

It is a simple, though tedious, calculation to verifyQ that this is indeed the spectral curve 
obtained by rotating the one given in ||. 

We have yet to examine the residue behaviour of the functions, and this will be of vital 
importance in what follows later. For t ~ B 

p{t-B) 2(t-B) (3.19) 

and 

1 4k 

h = r-^-^ = Us. (3.20) 

J 2w A(t — B) K ' 

The other relations then give fa ~ — 1/s, and hence 

T 3 ~rtdiag(l,-1,0) 

which identifies the representation formed by the residues as the irreducible one. As it will 
be needed later we list here the residue behaviour of all the functions at both ends of the 
interval, 



h ~ V2/s, h ~ V2/s, f 3 ~ 0/s, h ~ h ~ as s - (3.21) 

and defining the variable s = 2 — s then 

h ~ V2/s, / a ~ 0/s, / 3 ~ v^/s, f A ~ 1/s, / 5 ~ 0/s, as s - 2. (3.22) 



1 I thank Conor Houghton for checking this 
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4 Approximate twistor data 



In principle the equations were are concerned with for cyclic fc-monopoles are solvable in 
terms of abelian integrals of genus (k — 1) fll3"| . However, for k > 2, to explicitly extract 
from the general solution the one which satisfies all the required boundary conditions 
appears to be a highly non-trivial exercise. Furthermore, one of the main motivations for 
constructing the Nahm data is so that it can be used as input in the numerical ADHMN 
algorithm ||, to provide a visualization of the monopole energy densities. Even if the 
required explicit solution could be determined it is by no means clear that it would be in a 
form suitable for obtaining numerical values; there are at present no numerical algorithms 
available to compute the Riemann theta function for a surface of genus greater than one, 
unless it happens to be of a very special form which allows Weierstrass reduction theory to 
be applied [||]. With this in mind we construct, in this section, approximate twistor data 
(ie. Nahm data and spectral curves) for C3 symmetric 3-monopoles. 

One of the key points in applying the following method is that we know explicitly 
the twistor data for one member of the family ie. the tetrahedral 3-monopole. The one 
parameter family of monopoles we are searching for is a geodesic in the monopole moduli 
space and it is known |TIJ that the transformation between the monopole moduli space 
metric and the metric on Nahm data is an isometry. Since the Nahm data has poles it 
follows from these two facts that the residues at these poles must be constant with respect 
to the geodesic parameter. The upshot is that we know the residue behaviour explicitly 
for all members of the one-parameter family, it is given by ( |3.21| ) and fl3.22| ). 

As equations (|3.8|) are integrable they possess a particular solution which is a single- 
valued expansion around the s = pole ie 



(4.1) 



where the pole coefficients are given by (|3.21 ) as 

ai 1 = a/2, o 2 -1 = y/2, a 3 _i = 0, a 4 _i 



05,-1 = "I- 



(4.2) 



We now need to determine the number of arbitrary constants in the above particular 
solution. This can be done using Painleve analysis [[l], |14| as follows. We consider the two 
term expansion 

f { = «£± + b^ 1 (4.3) 
s 

for r an arbitrary integer, and linearise the equation (|3.8|) to obtain 



r 

r 

r -3 

y/2 

-V2 



-y/2 -2\/2 
2V2 y/2 



r-1 

r-1 









' " 




b 2 









h 


















b 5 








(4.4) 
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The Kowalevski exponents (they are not strictly resonances since 03,-1 = 0) are determined 
by the vanishing of the determinant A of the above matrix. We have that 



A = (r-3) 2 (r-2)(r + l)(r + 2) (4.5) 

so there are three arbitrary constants in our particular solution, corresponding to the 
three (counted with multiplicity) positive roots of A. Furthermore, we see that one of the 
arbitrary constants appears at the linear level in the s expansion and the remaining two 
at quadratic level. Explicit calculation reveals that these may be taken to be 01,1, ai 2 and 
03,2- We obtain that, to cubic order in s, 

ft = V2/ s + a M s + a 1)2 s 2 + -^f a i,i s3 

f 2 = V2/ s + a hl s - a h2 s 2 + -^- a i,i s3 
h = a 3t2 s 2 

1 2 

U= 1/s - V2ai,is + -j=a ly2 s 2 - -a 2 1A s 3 

1 2 

fc= + v^ai,^ + -^=a li2 s 2 + -a^s 3 . (4.6) 

This is going to be the form of the functions in our approximate Nahm data. Note that we 
terminate the expansion at 0(s 3 ), as this will turn out to be sufficient for our needs, but 
it is a simple matter to improve the accuracy of the approximate data by keeping higher 
order terms in the above. 

At this stage we have a candidate three-parameter family of approximate data, from 
which we need to select the correct one-parameter family. The constraining conditions 
arise from consideration of a symmetry of the equations ( |3.8| ). As before, let s = 2 — s, 
then if fi(s) is a solution we can construct a second solution gi(s) as 

[gi(s),g 2 {s),g 3 {s), gi {s),g 5 {s)} = [^(J), / 3 (s), / 2 (s), / 4 (s), -/ 4 (s) - / 5 (s)]. (4.7) 



Note that if fi(s) is a solution of the form ([4.6| ) with a pole at s = then gi(s) will have 
a pole at s = 2. Furthermore, the residues at this pole will be precisely those of ( |3.22p 
which we require. Thus, we take our approximate solution to be fi(s) for s G [0, 1] and 
gi(s) for s G [1,2]. There is now the matching condition at s = 1 that (?i(l) = /«(1), for 
i = 1,..,5. The first and fourth of these equations are identities, while the second and 
third are equivalent, so we are left with the two matching conditions 

/ 2 (1) = / 3 (1) and / 4 (1) = -2/5(1). (4.8) 

These give two relations between the three parameters 01,1,01,2,03,2 and thus determine 
the sought after one-parameter family. Using the series (|4.6|) the matching conditions give 
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that 

a 3 ,2 = V 2 + a 1;1 - a lj2 + ~l^- a i,i 

ai, 2 = ^{l-V2a 1A -la 2 11 ) (4.9) 
6 5 ' 

so that everything is determined by the parameter a 1A which we now relabel as a. In detail 
the functions are 

1_ ^--a^ + ^aV 
5 ; 20 

1_ ^--a^ + ^aV 
5 ; 20 

r- 2 2 7^/2 2 n 2 

- V2a a ) H a )s 

5 ; 20 ; 

- v^a - -a 2 )s 2 - -aV 
5 5 



/l = 


n 

V2/s + as + —( 

O 


f2 = 


V^/s + as- ^— ( 
o 


h = 


(V2 + a-^(l- 


h = 


l/s + V2as + -(l 


h = 


-1/s - V2as + ^ 



5 



(4.10) 



The next item to consider is the expression for the 'constants', a,/3 and 7. Of course 
they are only constant for exact solutions of the equations so for our approximate solutions 
( |4.10| ) they will have higher order corrections in s. Explicitly we find that a = + 0(s A ), 
= p Q + 0(s 2 ) and 7 = 7o + 0(s 2 ) where 



a = 12V2o 

1/ r- 29v/2 
A,= --(4^2 + 10a + — —a 2 ^ 



3 y 10 
70= |(5-5v / 2a-2a 2 ). (4.11) 

The second piece of approximate twistor data is thus the family of approximate spectral 
curves 

V 3 + a oV ( 2 + 7 oC 3 + ^o(C 6 + 1) = 0. (4.12) 

It is useful to analyse these curves, as it gives an indication as to the accuracy of the 
approximate construction. 

The first comparison that can be made is with the the exact tetrahedral 3-monopole 
spectral curve fl3.18|) . This has 

a = o, 7 = -5V2A _^ r(i/6)3r(i/ 3) 3 ^ 

For the approximate curve to have «o = requires a = 0, then by the above formulae we 
obtain 

«o = 0, 70 = -5^2$), -A) = 4V2/3 « 1.89. (4.14) 
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Note the surprising result that the relation between (3$ and 70 is exact in this case, so that 
the approximate spectral curve has tetrahedral symmetry. The value of (3q is also very close 
to the true j3 value, considering the low order approximation we chose. Clearly the accuracy 
of the approximate curves could be improved by performing an identical construction as 
above but taking a higher order approximating series. 

By considering asymptotic spectral curves we shall now determine the range of the 
parameter a, and also make some further comparisons. The asymptotic monopole config- 
urations along the exact geodesic can be determined from the corresponding rational map 
Q. At one end the configuration is that of three well-separated unit charge monopoles 
on the vertices of a large equilateral triangle in the xia^-plane. At the other end the con- 
figuration is asymptotic to an axisymmetric 2-monopole on the negative X3-axis, with the 
a^-axis the axis of symmetry, and a unit charge monopole on the positive X3-axis. These 
two clusters are well-separated with the distance from the origin of the 1-monopole being 
twice that of the 2-monopole. 

The spectral curve of a 1-monopole with position (xi,X2,x 3 ) is 

77 - (xi + ix 2 ) + 2x 3 ( + (x x - ix 2 )C = (4.15) 

and is called a star. By taking a product of such stars the asymptotic spectral curves 
corresponding to the above configurations can be determined ||. Taking three unit charge 
monopoles in the xiX2-plane with coordinates x\ + ix 2 = ibuj\ j = 0, 1,2 and to = e 27 ™' 3 
gives the product of stars 

= ( v -ib(l + ( 2 ))(ri-iku(l + ujC 2 ))(ri-ibu(u + C 2 )) = V 3 + vC 2 3b 2 - i(l + C 6 )& 3 . (4.16) 

Hence the asymptotic curve has dihedral D 3 symmetry, since 7 = 0, with the other two 
parameters, a positive, (3 negative, satisfying the relation 

a(-(3)- 2 / 3 = 3. (4.17) 

We now define the upper limit a + of the parameter a to be that for which the approximate 



curve has D$, symmetry ie. 70 = 0. By equation (4.11) this determines a + to be 



a+ =^^«0.60. (4.18) 

At this value of a it can be calculated that 

ao(-A)r 2/3 « 3.8 (4.19) 

which should be compared with (|4.17|) . 

The product of stars of a 1-monopole at position (0, 0, h) and an axisymmetric 2- 
monopole at (0, 0, —6/2) is 

= (, + - ^ + + £ K>) = f + - **) + ^ + £ )C* (4.20) 
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Hence the asymptotic curve is axisymmetric about the X3-axis, since j3 = 0. For large 
separation the asymptotic relation is that a is negative and 7 is positive with 



«7" 2/3 = « -1.9. (4.21) 

In a similar fashion to above, the lower limit a_ is defined to be the value of a for which 
the approximate curve has axial symmetry ie. (3q = 0. This gives 

a_ = ^(3^ - 25) « -0.89 (4.22) 
29 

and at this parameter value 

-2/3 _ 

(15(65 - 2V/5)) 2 / 3 



= „J (87 ''1,„. - -1.7 (4.23) 



which is to be compared with ( |4.21| ). 

In Figure 1, we plot the approximate spectral curve coefficients, a >Ab7o f° r a £ [ a -> 
Having discussed the approximate spectral curves it is now time to use the approximate 
Nahm data for its intended purpose. We take it as input for the numerical ADHMN 
construction developed in an earlier paper ||. Figure 2 shows the output, in the form of a 
surface of constant energy density, for each value of the input parameter a. The surfaces 
shown correspond to the values a = 0.4, 0.2, 0.1, 0.0, —0.1, —0.2, —0.4. The low energy 3- 
monopole scattering process is as follows. Initially, Fig 2.1, there are three unit charge 
monopoles on the vertices of a contracting equilateral triangle in the xiX2-plane. As they 
merge each monopole raises an arm, so that they link in a kind of maypole dance, Fig 2.2. 
The legs of the monopoles continue towards the centre, Fig 2.3, until they too merge and 
the tetrahedral 3-monopole is formed, Fig 2.4. Next the top segment of the tetrahedron 
separates from the bottom, Fig 2.5, and as it moves up the 23-axis it leaves behind a torus 
with three prongs, Fig 2.6. As the single monopole continues its journey up the X3-axis it 
becomes more spherical, and the 2-monopole smooths out into a torus as it moves down 
the X3 axis. 



5 Higgs zeros 

Numerical evidence suggests Q, 12 1 that the tetrahedral 3-monopole has five zeros of the 



Higgs field, despite the fact that it is a charge three monopole. A scattering geodesic 
through the tetrahedral 3-monopole has been investigated in detail and the associated 
dynamics, creation and annihilation of the Higgs zeros tracked numerically 0. Since the 
cyclic scattering of three monopoles passes through the tetrahedral 3-monopole, we have an 
opportunity to study further this novel phenomenon of extra Higgs zeros. We compute an 
analytical approximation to the Higgs field, and find that it is consistent with a conjectured 
behaviour of the Higgs zeros. 
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(5.3) 



Finding the Nahm data effectively solves the nonlinear part of the monopole construc- 
tion but in order to calculate the Higgs field the linear part of the ADHMN construction 
must also be implemented. Given Nahm data (Ti,T2,T 3 ) for a fc-monopole we must solve 
the ordinary differential equation 

(l2fc^ + Ik ® XjVj + iTj ® (Tj)v = (5.1) 

for the complex 2k- vector v(s), where 1*. denotes the k x k identity matrix, <jj are the 
Pauli matrices and x = (xi,x 2 ,x 3 ) is the point in space at which the Higgs field is to be 
calculated. Introducing the inner product 

(vi, v 2 ) = / vjv 2 ds (5.2) 
jo 

then the solutions of ( |5.1|) which we require are those which are normalizable with respect 
to (|5.2|) . It can be shown that the space of normalizable solutions to (|5.1|) has (complex) 
dimension 2. If v 1; v 2 is an orthonormal basis for this space then the Higgs field $ is given 
by 

((s-l)vi,vi) ((s-l)vi,v 2 ) 
((s-l)v 2 ,vi) ((s-l)v 2 ,v 2 ) 

The strategy adopted is to use the approximate Nahm data of the previous section and 
compute approximate solutions of (|5.1|) . 

To begin we consider the initial value problem of fl5.1|) at the pole s — 0, which has the 
form 

= B s v (5.4) 

as 

where B s is a regular 6x6 matrix function of s G [0, 2). This is a regular-singular problem 
and the eigenvalues of B are {1, 1, 1, 1, —2, —2}, hence there is a four-parameter family of 
solutions which are regular for s G [0, 2). Expressing the solution as a series 

oo 

v(s) = ^d^V (5.5) 

3=1 

the four arbitrary parameters may be taken to be d\, .., where 

dW = (di, da, d 3) d 4 , iV2d 2 , -iV2d 3 f. (5.6) 

Similarly for the intial value problem at s = 2 there exists a four-parameter family of 
solutions v(s) which are regular for s G (0, 2] and which may be expressed as a series in 
8 = 2 - s 

oo 

v(s) = 5^h^'V. (5.7) 

In this case the four arbitrary parameters are contained in 

h (1) = (h h h 2 , h 3 , iV2h h hi, th 3 /V2Y. (5.8) 
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Taking the solution to be v(s) for s G [0, 1] and v(s) for s G [1, 2] results in the matching 
condition v(l) = v(l) at s = 1. Since this is a vector equation with six components this 
gives six constraints involving the eight arbitrary parameters. We thus obtain the required 
two parameter family of solutions from which to construct an orthonormal basis. 

To perform the approximate ADHMN construction we use the approximate Nahm data 
(with parameter a) of the previous section and truncate the series ( |5.5|) and (|5.7|) at quartic 



order. Since all functions are now simple series the required integrals are elementary. The 
above procedure was implemented using MAPLE. Despite the low order approximations 
used the calculations are very involved, and to simplify matters we restrict to calculating 
the Higgs field on the x 3 -axis only. That is, we set x\ = x% = and £3 = L in the above. In 
this case a convenient choice of gauge (equivalent to a choice of orthonormal basis) exists 
in which the Higgs field is diagonal 



<f 
-<p 



(5.9) 



The output of the MAPLE program then consists of the function (p, which depends on 
the coordinate along the 23-axis L and the input parameter a, which determines which 
monopole configuration along the geodesic we are studying. Even in this simplified case 
the result is rather complicated (the expression fills a full page!). We plot, in Figure 3, p 
as a function of L for four values of a. The first (Figure 3.1) is for a = 0, which is the 
tetrahedral 3-monopole. It can be clearly seen that there are two distinct zeros of the Higgs 
field along this line. One is at the origin and the other is at L « 1, which is associated with 
one of the vertices of the tetrahedron. This result is in good agreement with the numerical 
results of [0, |i~2f . It confirms that the tetrahedron has four zeros, each with a local winding 
+1, on the vertices of a tetrahedron and a fifth zero, which we refer to as an anti-zero since 
it has local winding —1, at the origin. 

The remaining plots in Figure 3 are evidence to support the following conjecture on 
the dynamics of the Higgs zeros. Initially the three monopoles are well-separated, so there 
can only be three zeros. They are positioned on the vertices of an equilateral triangle in 
the Xix 2 -plane. At all times the zeros of the Higgs field must be consistent with the C 3 
symmetry, which leads us to the following description. 

The three zeros on the vertices of a triangle approach the a^-axis, but fall below the 
origin. During their approach there is a critical point on the o^-axis where there is a zero 
anti-zero creation event, sending an anti-zero down the X3-axis and a plus zero up the 
X3-axis (Figure 3.2, a = 0.05). The anti-zero reaches the origin (this is the tetrahedral 
3-monopole, Figure 3.1, a = 0) and continues down the x 3 -axis (Figure 3.3, a = —0.05) 
until it eventually meets up with the three positive zeros as they hit the X3-axis. Then we 
have lost all anti-zeros and are left with the +1 zero of the unit charge monopole moving 
up the X3-axis and the +2 zero of the axisymmetric 2-monopole moving down the a^-axis. 
As a last check we verify that for a sufficiently large there are no zeros on the o^-axis 
(Figure 3.4, a = 0.2). 

From Figure 3 we see that the number of zeros on the a^-axis has a simple interpretation. 
The shape of the ip curve remains the same, but as a decreases the curve moves up the 
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plane. Thus the number of roots of this curve increases from zero to an instantaneous 
double root and finally to two simple roots, with the separation between the two roots 
then monotonically increasing. 

We now address the question of the uniqueness of the above interpretation of the motion 
of the Higgs zeros drawn from the data. It turns out that there is a second possibility which 
is more complicated than the first and involves a splitting into anti-zeros of each of the 
three initial Higgs zeros. The results given so far are consistent with both conjectures and 
an additional check must be made to decide between the two. For the first possibility, 
which will turn out to be the correct one, we have seen that for configurations just prior 
to the formation of the tetrahedral 3-monopole there are two zeros on the X3-axis which 
have a local winding number +1 and —1 respectively. An analysis of the second possibility 
demands that these zeros have local winding numbers +1 and +2. To settle the issue a 
numerical calculation is performed. First we select a suitable configuration, which is taken 
to be a = 0.02. Figure 4 shows a plot of the length squared of the Higgs field (solid line) 
and a component of the Higgs field (dashed line) along the axis x% = L. From this the 
positions of the two zeros can be approximately read off as L « 0.6 and L « 0.15. Next we 
apply the numerical algorithm given in to compute the local winding number Qr(L) of 
the normalized Higgs field on a two-sphere of radius R, with centre (xi, x 2 ,x 3 ) = (0, 0, L). 
As a first check we compute that £^2(0) = +3, so that when all zeros are inside the chosen 
two-sphere the local winding number indeed gives the number of monopoles. Next we 
examine the zero which is highest on the X3-axis and confirm that <3q.i(0.6) = +1. Finally, 
the computation that (5o.i(0.15) = — 1 shows that the second zero on the X3-axis is an 
anti-zero and thus rules out the second possibility for the motion of the Higgs zeros to 
which we alluded above. 

For geodesies involving monopoles which are linearized on a quotient curve which is 
elliptic, it has been conjectured |7|, that the special 'splitting points', where anti-zeros 
appear/disappear, are associated to singular points where an elliptic curve becomes ratio- 
nal. Presumably when the geodesic involves monopoles which are linearized on a higher 
genus surface (as here) the 'splitting points' will correspond to vanishing cycle points of 
the surface. 

6 Cyclic 4-monopoles 

For C4 symmetric 4-monopoles there are three types of geodesies, corresponding to the 
surfaces £4, £4 and £|. The first of these is the essentially planar scattering where the 
four monopoles instantaneously form an axially symmetric configuration and subsequently 
scatter through an angle of 45°. We shall not be concerned with this type of scattering. 
The other two processes are fully three-dimensional and are of more interest. The £4 
geodesic has an increased symmetry from C4 to -D4 and describes a scattering through the 
cubic 4-monopole [[J. It is possible to investigate this geodesic in exactly the same manner 
as was done earlier for the C3 geodesic through the tetrahedral 3-monopole. However, the 
£4 geodesic does not contain any Platonic monopole configurations, and so at first sight it 
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appears that we require a new ingredient to deal with this case; since we do not know the 
exact Nahm data at any point on the geodesic we can not read-off the residue behaviour. 
It turns out that it is instructive not to deal with each case separately, but to apply a 
unified treatment. This is presented in the following. 

First of all, we use equation ( |3.4| ) to obtain the form of the Nahm data from the algebra 
It is convenient to change notation and write 



Ti = - 
2 



" 


-Qi 





9o 


Ql 













Q2 







. -Qo 





?3 






To 



qi g 

<?i q 2 

qi g 3 

<?o q 3 



pi 

p 2 

p 3 

po 



(6.1) 



Nahm's equation then becomes the set of equations (equivalent to the 4-particle Toda 
chain) 

Pi = 2 (?? - 9?-i)> <?j = 9i(Pj - Py+i) (6.2) 

where the indices are to be read modulo 4, in accordance with the periodicity of the chain. 
Since we are considering strongly centred monopoles we also impose the standard Toda 
chain zero momentum relation 



5>i = °- 

3=0 



(6-3) 



We shall generally use this relation to eliminate po from the equations, but for certain 
purposes it is useful to preserve the symmetry of the system and not explicitly solve for pq. 
The spectral curve of the Nahm data ( |6.1| ) is 



rf + arfC 2 + P(( 8 + 1) + l( 4 + = 
where the real constants are 



(6.4) 



a = + 2 Pj( 2 Pi+i + P3+2)} 
3=0 

P = - fin 

j=0 

3 ^3 

7 = ^I[Pj + ^E{^(<?!+2 + 8p j - 1 p j+2 )} 

j=0 z j=0 

3 

6 = £{9?(Pi-i +P3-2) + 8 PjPj+iPj-i}- 

3=0 



(6.5) 
(6.6) 
(6.7) 
(6.8) 
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We now examine the residue behaviour at s = 0. Putting 



1j ~ r j/ s i Pj ~ -Rj/s as s 



(6.9) 



into equation 



gives 



2R 3 + r) - r)_ x 

Ra 




0. 



r^Rj-Rj+x + 1) 
Substituting the first of these equations into the second we obtain 

^(^ 2 +i-2r| + r|„ 1 + 2) 



0. 



(6.10) 
(6.11) 



(6.12) 



Now if all the rj are non-zero then, recognising the bracket as the second order difference 
approximation to the second derivative, the general solution of the difference equation is 



r j = c + cij - 2j" 



(6.13) 



for constants c ,Ci. However, this solution does not satisfy the periodic boundary condi- 



and thus we have proved that at least one of the r,'s must be zero. 



tions rj + 4 

By symmetry we can choose r = 0. All the R^s must be non-zero (in order for 
the matrix residues to form the irreducible 4-dimensional representation of su(2)) which 
together with equation ( 6.10|) implies that r\r% ^ 0. Furthermore, r 2 7^ 0, since if r 2 = 
then the solution is r\ = r\ = 1 which by equation (|6.10| ) determines (i? , Ri, R2, R3) = 
(1/2, —1/2, 1/2, —1/2) giving the reducible representation 2© 2. Now that we have r = 
and rir 2 r 3 7^ the equations simplify to the linear system 



(6.14) 



Note that the matrix occuring in the above is the Cartan matrix of the Lie algebra A3. 
The unique solution is 



2 -1 







r\ 




' 2 " 


-1 2 


-1 




rj 




2 


-1 


2 




. r l . 




2 



r = 0, 



4. 



which gives from equation fl6.10j ) that 

(Ro, R u R 2 , R 3 ) = (3/2, -3/2, -1/2, 1/2) 



(6.15) 



(6.16) 



which identifies the matrix residues as forming the required irreducible representation 4. 

The upshot is that the only freedom in the choice of residues is which one of the r^'s 
is chosen to be zero. There are some choices of sign available in taking the square roots 
in ( |6.15[ ) but these only lead to an inversion of the whole configuration and are of no 
importance. By a suitable choice of gauge it is always possible to choose r = and so at 
first glance it may appear a puzzle as to how the different types of geodesic are selected. 
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However, there is a similar second pole at s = 2 at which again one of the associated 
residues fj must be zero. Since we are not using a basis in which the symmetry property 
Tj(2 — s) — Ti(s) 1 of the Nahm data is manifest, then once the gauge freedom has been 
used to fix ro = then we have no freedom left to specify which rj vanishes. The different 
geodesies are distinguished by having different r/s being zero. Another way to view this 
is that the explicit symmetry transformation relating Nahm data at s to Nahm data at 
2 — s will be different for different geodesies in the same basis. Of course, by general 
arguments it is true that a basis exists in which the symmetry property is the manifest one 
given above, but the important point is that this basis is not the same one for different 
geodesies. This means that there is no point looking for a general "good basis" at the start 
of the calculation, since this will vary with each solution. It does however mean that all 
the geodesies can be treated at once and the different cases identified near the end of the 
calculation by the imposition of the different symmetry properties. 

We now construct the approximate twistor data, using the method introduced earlier. 
A particular solution, which is a single- valued expansion about the s = pole with residues 
Tj, Rj as determined above, is constructed for the functions qj,Pj- We omit the details since 
the method is just as in the C3 case considered earlier. Applying Painleve analysis gives a 
linear system with the matrix determinant 

A = 8(r + l)(r + 2)(r + 3)(r - 2)(r - 3)(r - 4) 2 (6.17) 

giving the Kowalevski exponents. The solution thus contains four arbitrary constants, 
which is the required number to coincide with the constants a, /3, 7, 5 in the spectral curve. 

Again we work with expansions to cubic order in s and take these as the approximate 
Nahm data for s G [0,1]. To determine a one-parameter family of monopoles we must 
impose a symmetry relation to construct the Nahm data for the second half of the interval 
ie. s e [1,2]. 

The case E4 

For this geodesic the symmetry property of equations (|6.2|) which we exploit is the 
transformation under s 1— > 2 — s of 

(90,91,92,93) ^ (92,9i,9o,93) (6.18) 
(Po,Pi,P2,P3) >-> (pa,P2,Pi,Po)- (6.19) 

Using this transformed data in the second half of the s interval gives the matching condi- 
tions at s = 1 

9o (l) = 92(1), Pi(l) = P2U), Pi(l) + Pa(l) + 2p 3 (l) = 0. (6.20) 

Applying these three constraints to the four-parameter family of solutions we arrive at 
one of the sought after one-parameter families. The constant terms in the spectral curve 
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coefficients are given by 

a = 20V3a (6.21) 
15 

Pa = — + 7V3a + 4a 2 (6.22) 

7o = 105 - 7(h/3a + 80a 2 (6.23) 
5 = 0. (6.24) 

In fact it is easy to see that the imposed symmetry forces the reductions q\ = and 
£>3 = — p2 which in turn gives 5 = (not just 5q = 0). Thus the approximate data has 
exact -D4 symmetry. We therefore already know that we are considering either the £4 or 
the £4 geodesic. The confirmation that we have the second of these two possibilities is 
provided by examination of the approximate spectral curve for a = 0. This curve is given 
by 

« = 0, 7o = 14A,, A) = 15/2. (6.25) 

The first two of these relations are enough to ensure that this curve has cubic symmetry. 
Thus we have a geodesic containing a cubic configuration and hence are considering the 
£f geodesic. For comparison the exact cubic 4-monopole has a curve given byQ 

a = 0, 7 = 14/3, f3 = 3r(l/4) 8 /(10247r 2 ). (6.26) 

Because of the way that (3 enters into the spectral curve it is really the fourth root of (3 
which determines the length scale of the cubic monopole. Thus the appropriate comparison 
is 

/3 1/4 = 1.73, /? 1/4 = 1.65 (6.27) 

which again demonstrates the reasonable accuracy of such a low order approximation 
scheme. 

Using the same kind of analysis as in the C3 case considered earlier, it is possible to 
determine the range of the parameter a by examining the limiting spectral curves. We 
apply the numerical ADHMN construction to this approximate Nahm data for the values 
a = 0.4,0.2,0.1,0.0, —0.1, —0.2, —0.4. The resulting energy density surfaces are shown in 
Figure 5. The four monopoles approach on the vertices of a contracting square Fig 5.1, 
and bend as they merge Fig 5.2, until they form the cube Fig 5.4. The top and bottom 
portions of the cube then pull apart Fig 5.5, and become increasingly toroidal as the pair 
of charge 2 monopoles separate along the a^-axis Fig 5.7. 

The case £4 

For this geodesic the relevant symmetry property of equation ( |6.2|) is the transformation 
under s 1— > 2 — s 

(go, qi, 92, qa) | -> (lu Qo, qa, Qs) (6-28) 

(P0,Pl,P2,P3) >-> (P2,Pl,P0,P3) (6-29) 



2 There is a factor of 16 error in ref. || 
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leading to the three matching conditions 

q (l) = gi (l), q 2 (l) = q 3 (l), Pl (l) + 2p 2 (l) +p 3 (l) = 0. (6.30) 
In this case the approximate spectral curve coefficients are a little more complicated 
a = 20V3a (6.31) 
30v / 3 + a(138 + 72 v / 3) + a 2 (36 + 25 v / 3) , c nn . 

Po = (6 - 32) 

-1260 - 1680\/3 + a(3360 + 2100^) + a 2 (4530 + 2500^) 
70 ~ (v / 3 + 2)(9 + 5 v / 3) ' 

-2400 - 1200\/3 + a(2400 + 1360x/3) + a 2 (520 + 320^) , M) 

° ~ (v / 3 + 2)(9 + 5 v / 3) ' ' 

The salient feature is that 5q ^ 0, so the symmetry is C4 and not D4, thereby identifying 
the geodesic as the one associated with the surface £4. 

In figure 6 we show energy density plots for the parameter values 
a = 0.6, 0.5, 0.4, 0.3, 0.2, 0.0. This scattering appears very similar to the C 3 scattering 
discussed earlier. The four monopoles approach on the vertices of a contracting square Fig 
6.1, and as they merge they link arms forming a pyramid shaped configuration Fig 6.3. 
The top of the pyramid breaks off, Fig 6.5, and travels up the a;3-axis, while the charge 
three base travels down the X3-axis. As they continue to separate, Fig 6.6, the unit charge 
monopole becomes more spherical and the base deforms closer towards the axisymmetric 
3-monopole. 



Numerical evidence suggests [12| that the cubic 4-monopole does not possess anti- 
zeros, in contrast to the tetrahedral 3-monopole. Figures 5 and 6 appear to add some 
understanding to this result, since it is the £4 geodesic, rather than the £ 2 geodesic, which 
appears to be the closest 4-monopole analogue of the £3 geodesic. Hence we might expect 
the pyramid-like 4-monopole, rather than the cubic 4-monopole to be the one which has 
anti-zeros. 



7 Conclusion 

We have proposed an analytical method to obtain approximate Nahm data which, when 
combined with the previously introduced numerical ADHMN construction, provides an 
efficient approximation scheme for the construction of monopoles. This has been applied 
to the study of monopoles with cyclic symmetry and shown to produce good results. The 
scattering processes reveal exotic dynamics, and indicate that the key to understanding 
this type of scattering lies with a more detailed knowledge of the zeros of the Higgs field. 
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Figure Captions 

Fig 1. Approximate spectral curve coefficients 

Fig 2. C3 symmetric energy density surfaces. The parameter values are a = 

0.4,0.2,0.1,0.0,-0.1,-0.2,-0.4 

Fig 3. Plot of the Higgs component ip on the rr 3 -axis, for configurations with (l)a = 
0; (2)a = 0.05; (3)a = -0.05; (4)a = 0.2. 

Fig 4. Plots of the length squared of the Higgs field (solid line) and a Higgs component 
(dashed line) along the x 3 -axis for a = 0.02 

Fig 5. D4 symmetric energy density surfaces. The parameter values are a = 
0.4,0.2,0.1,0.0,-0.1,-0.2,-0.4 

Fig 6. C 4 symmetric energy density surfaces. The parameter values are a = 
0.6,0.5,0.4,0.3,0.2,0.0 

Note: Figures 2,5 & 6 are attached gif files cyclicfig2.gif etc. 
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Figure 2: View figure by typing xv cyclicfig2.gif 
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